Localized states in sheared electroconvection 
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Abstract. - Electroconvection in a thin, sheared fluid film displays a rich sequence of bifurcations 
between different flow states as the driving voltage is increased. We present a numerical study 
of an annular film in which a radial potential difference acts on induced surface charges to drive 
convection. The film is also sheared by independently rotating the inner edge of the annulus. This 
simulation models laboratory experiments on electroconvection in sheared smectic liquid crystal 
films. The applied shear competes with the electrical forces, resulting in oscillatory and strongly 
subcritical bifurcations between localized vortex states close to onset. At higher forcing, the flow 
becomes chaotic via a Ruelle-Takens-Newhouse scenario. The simulation allows flow visualization 
not available in the physical experiments, and sheds light on previously observed transitions in 
the current- volt age characteristics of electroconvecting smectic films. 
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Driven, dissipative nonlinear systems sometimes exhibit 
spatially localized structures which are analogous to soli- 
tons [1]. Examples are found in systems as diverse as veg- 
etation patterns [2], vibrated granular media [3] and fer- 
rofluids [4]. Such states can also arise in fluid- mechanical 
settings such as binary fluid convection [5,6] (where they 
have been called "convectons" [7]), in electroconvecting 
nematic liquid crystals [8] , and in very general model equa- 
tions [9,10]. In this Letter, we describe a new and unex- 
pected type of localized patterns: solitary vortex states 
in two-dimensional, sheared electroconvection. We show 
numerically that these states can be extremely localized, 
consisting of only a single, isolated vortex surrounded by 
a uniform background state. The presence of localized 
states in this system is surprizing and interesting because 
it results from the interaction of a circular Couette shear 
with two-dimensional convection in an especially simple, 
highly symmetric geometry. 

Theses states can be experimentally realized using thin 
free-standing films of smectic liquid crystals [11-16]. Here, 
direct numerical simulation allows us to study the spa- 
tial structure of the full velocity, charge and potential 
fields [16]. This approach compliments existing theory [12, 
16] and experiments [11-15] on this system, which mainly 
consisted in observations of the total current through the 
thin film, without flow visualization. The simulation re- 
veals localized states in the form of vortices which travel in 



the direction of the applied shear. These are preceded by 
lower- amplitude, extended traveling and oscillatory vortex 
states. At sufficiently high levels of electrical forcing, the 
flow becomes chaotic via a Ruelle-Takens-Newhouse sce- 
nario [17]. Our results serve to strongly motivate new ex- 
periments and theory to elucidate the dynamics of sheared 
2D convection. More generally, this system with naturally 
periodic boundary conditions, and for which forcing and 
shear are independently controllable, will be an interest- 
ing place to examine recent ideas on higher dimensional 
invariant manifolds on the way to turbulence [18,19]. 

Our numerical study simulates a laboratory experiment 
shown schematically in Fig. 1. The system consists of a 
submicron-thick liquid crystal film freely suspended be- 
tween concentric circular electrodes. The weakly conduct- 
ing film is driven to convect when a sufficiently large elec- 
tric potential is imposed across it. The inner edge of the 
annular film is held at absolute potential V with respect to 
infinity, while the outer edge is held at zero potential. In 
addition to the control parameter associated with this elec- 
trical driving force, it is possible to independently rotate 
the inner electrode, which imposes an azimuthal Couette 
shear on the film. The experimental signature of convec- 
tion consists of measurements of the total current through 
the film, which is increased by convective flow. 

The film develops a surface charge configuration which 
is unstable to the applied voltage. This instability is 
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Fig. 1: (Color online) Schematic of the sheared 2D annular 
film. 

closely analogous to that of Rayleigh-Benard convection, 
in which an inverted density distribution is unstable 
to buoyancy forces. As in Rayleigh-Benard convection, 
there are two important dimensionless parameters. The 
Rayleigh-like number 71 describes the ratio of the electri- 
cal forcing to the viscous and electrical dissipation; this 
serves as the main control parameter. The Prandtl-like 
number V is the ratio of charge relaxation time r q to vis- 
cous relaxation time r v . The annular geometry of the film 
is characterized by the radius ratio a. These parameters 
are given by [12] 

elV 2 e r] n 

71 = — — , V = and a = — , (1) 

ar] pad r Q 

where p, 77 and a are the 2D mass density, shear viscosity 
and electrical conductivity, respectively. The width of the 
film is d = r Q — r^, and eo is the permittivity of free space. 

Rotation of the inner electrode introduces a circular 
Couette shear to the base state below the onset of con- 
vection. The strength of the shear is characterized by a 
dimensionless shear Reynolds number 

Tie = pdouri/rj, (2) 

where uj is the angular rotation frequency of the inner elec- 
trode. The basic equations, given in detail in Refs. [12] 
and [16], are nondimensionalized using the film width 
d as the unit of length and the charge relaxation time 
r q = eod/ a as the unit of time. Recall that a is the 2D 
conductivity. In these units, the dimensionless rotation 
rate of the inner electrode is Q = cor q . 

The dimensionless electrical current through the film is 
described by the Nusselt number, Nu = I/Io, where Iq is 
the fraction of the current contributed by pure conduction 
at the same voltage. In a convecting film, Nu > 1. The 
quiescent film undergoes a bifurcation from conduction to 
convection when a voltage larger than a critical value V c 
is applied, corresponding to a critical value 71 C1 which in 
general depends on ft, P and Tie. Typically, Tl c ~ 100. 
Just above onset, the flow is characterized by a number 
of counter-rotating vortices arranged around the annulus, 
filling its entire circumference. Under shear, the vortex 
pattern rotates, traveling in the same sense as the inner 
electrode. 



We generalize our previous numerical scheme [16] to en- 
compass the case of nonzero shear. Denoting the quanti- 
ties in the base state below the onset of convection by 
a superscript zero, the only new nonzero quantity that 
appears in the perturbed equations (Eqns. 18 and 19 in 
Ref. [16]) is the radial derivative of the base state stream 
function <j)^ given by 

dr 1 - a 2 V r(l -a) 2 )' K J 

where r is the dimensionless radial coordinate in units of 
d. With this change, the numerical treatment described 
in Ref. [16] can be carried out for the sheared case. 

The direct numerical simulations used a pseudo-spectral 
scheme, based on the Fourier Galerkin method in the az- 
imuthal direction and the Chebyshev collocation method 
in the radial direction [20]. The governing equations, 
given in full in Ref. [16], are the 2D Navier-Stokes equa- 
tion, the continuity equations for mass and charge, and 
one Maxwell equation for the nonlocal, coupled relation- 
ship between the charge and electric potential on the film. 
We used the streamfunction and the vorticity as primitive 
variables. These were computed in the annular geometry 
< < 2tt and a/(l - a) < r < 1/(1 - a). The bound- 
ary conditions are periodic in the direction, while a no- 
slip boundary condition holds at the edges of the annulus. 
To make sure the computations are well resolved, we em- 
ployed up to 32 Fourier modes, up to 40 radial Chebyshev 
collocation points and a variety of time- step sizes. We ver- 
ified the simulation results using different expansion orders 
and different time-step sizes [16]. 

Unlike in the case of 3D Taylor- Couette flow, here the 
2D applied shear is itself stable and always stabilizes the 
conduction state, suppressing the onset of electroconvec- 
tion. At the primary bifurcation at the onset of convec- 
tion, the critical azimuthal Fourier mode with mode num- 
ber m c becomes unstable. We denote lZ c (JZe = 0) = 71® 
and m c (Tle — 0) = m®. The simulation results for 71® 
and m® show excellent agreement with previous linear 
stability analyses [12,16]. Under applied shear, we find 
7l c > 71® and m c < m®. This suppression of convec- 
tion is described by a reduced critical control parameter, 
e = (7l c /7l®) — 1. This reduced quantity can be extracted 
from experimental current data with essentially no ad- 
justable parameters, and calculated using linear stability 
analysis [12]. The numerical simulation also produces val- 
ues of the reduced Nusselt number n = Nu — 1 that can be 
directly compared to experimental data in the nonlinear 
regime. We discuss such comparisons below. 

In addition to the integrated quantities like Nu, the 
simulation also provides directly the full vorticity, stream 
function, surface charge, and electric potential fields. 
These basic fields are not directly accessible to experi- 
ment, except by using rather invasive flow visualization 
techniques which tend to unduly perturb the system [11]. 
Thus, we can use the simulation to provide important new 
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insights into the flow dynamics as 1Z and IZe vary, and 
our results shed new light on some previously unresolved 
experimental observations, while also making interesting 
new predictions. 

The simulations clearly show that the vicinity of the on- 
set of convection under shear is much more complex than 
previously suspected, with new low-amplitude and oscilla- 
tory states predicted. The simulations show that all pri- 
mary bifurcations are continuous, or "supercritical" , even 
under shear. This is in contrast with previous experimen- 
tal studies [13, 14] which were interpreted to show that 
the bifurcation became hysteretic and discontinuous, i.e. 
"subcritical" , under a sufficiently large shear. This dis- 
crepancy may be traced to the fact that the excess current 
due to convection just above onset is very small, and can 
just barely be distinguished from the background conduc- 
tion current due to noise and drift effects. 

Fig. 2a shows raw experimental current-voltage data 
from Fig. 2 in Ref. [14] and the corresponding measured 
dimensionless quantities Nu — 1 vs. 1Z. Fig. 2c shows the 
results of a full numerical simulation at the same dimen- 
sionless parameters as the previous figure. The integrated 
quantity Nu — 1 found numerically is in reasonable agree- 
ment with the experimental data. However, an examina- 
tion of the flow pattern reveals several new features which 
cannot be deduced from the experimental data alone. 

Convection begins at a lower voltage than previously 
thought, and the first state encountered has a very low am- 
plitude, corresponding to a small value of Nu — 1. In the 
analysis of experimental data like that shown in Fig. 2a, 
the critical voltage V c was misidentified to be near the 
jump at Vj instead of somewhere close to V* . The nu- 
merical simulation shows that there is a sudden large in- 
crease in the current near Vj which is associated with a 
sharp change in the mode structure, while the true onset 
of convection is lower, near V* . Just a trace of the low am- 
plitude signal may be seen in the measured Nu — 1 shown 
in Fig. 2b. In some unpublished data [21], the effect was 
more pronounced than in Fig. 2b, and was suspected to 
be due to low amplitude supercritical states, though the 
matter could not be experimentally resolved at that time. 

The misidentification of V c means that the measured 
suppression of the onset on convection, represented by the 
experimental value of e, must be systematically too large. 
The difference in V* and V/ is typically of the order of the 
width of the hysteresis loop which is about 10 — 30% of V c , 
depending on other parameters. The systematic overesti- 
mate of the experimental suppression may be seen in Fig. 7 
of Ref. [13], where in the experimental data lie above the 
theoretical suppression from linear stability calculations. 

As 1Z gradually increases beyond the steady sheared 
convection regime, the simulation shows that a Hopf bi- 
furcation sets in leading to a periodic motion in a small JZ 
range, as shown in Fig. 2c. In addition to the usual rota- 
tion of the vortex pattern, the amplitude becomes oscilla- 
tory in time, while generally still filling the circumference. 
The steady sheared convection state does not show a time 
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Fig. 2: (Color online) (a) Experimental raw current- voltage 
data from Ref. [14] for a = 0.56, V = 76, and IZe = 0.124. 
The open (filled) triangles are for increasing (decreasing) volt- 
age, (b) The corresponding measured dimensionless connec- 
tive current, Nu — 1, vs. dimensionless control parameter 
e = (V/Vc) 2 — 1, extracted from the data in (a), (c) Sim- 
ulation results for Nu — 1 for the same parameters as in the 
previous experimental figure. The film undergoes a sequence of 
bifurcations as 1Z is slowly increased and then decreased (y). 
The sequence of flow regimes are conduction (*), steady ro- 
tating convection (A), and oscillatory rotating convection (o), 
with the vertical bars showing the amplitude of oscillation. A 
large hysteretic bifurcation is located at the jump in the di- 
mensionless current that signals changes in the structure of 
the azimuthal Fourier modes. 

dependence in Nu — 1 until these amplitude oscillations 
set in. 

Oscillatory motion under a moderate shear correspond- 
ing to IZe = 0.124 begins at about 1.18 x 1Z C . Fig. 3a 
shows numerical time-series data of the oscillating con- 
vective current at 1Z = 371.3. After ~ Wr q , we find a 
one-frequency limit cycle is established. The correspond- 
ing phase space reconstruction of this cycle is shown on 
the right of Fig. 3a, using the time-delay method with a 
sufficiently long time series of Nu — 1. The oscillatory 
regime spans the range 365 < 1Z < 388. Fig. 3b shows 
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Fig. 3: (Color online) Numerical data for the oscillatory con- 
vective current Nu — 1 for a = 0.56, V — 75.8, and IZe — 0.124: 
(a) the time-series data for 1Z = 371.3, together with the cor- 
responding phase space trajectory, reconstructed using a time- 
delay method; (b) the fundamental frequencies (•) and the 
oscillation amplitudes (■) of Nu — 1 as a function of 1Z. 

the oscillation frequencies and amplitudes over this range. 
The oscillation frequency for 1Z = 371.3 is 0.40 in units of 
r~ x . Simultaneously, the rate of rotation of the vortex pat- 
tern itself, which is carried around by the imposed shear, 
measured in the range 1Z = 390 — 394, was 2.36 rad/r g 
counterclockwise, corresponding to a physical frequency of 
0.38 Tq 1 . Thus, the frequencies of the amplitude oscilla- 
tion and of the overall pattern rotation are comparable in 
the narrow window of oscillatory convection. 

The oscillatory convection regime is not always present 
under shear, however. We did not observe oscillations for 
very small shears, such as IZe = 0.01 and IZe = 0.06, for 
similar supercritical values of 1Z. 

Previous laboratory experiments [13, 14] measured the 
time averaged total current passing through the film, pro- 
ducing data like that shown in Fig. 2a. A re-examination 
of the time series that were averaged over, close to the 
sudden current increase that signals the end of the oscil- 
latory regime, failed to unambiguously identify an oscilla- 
tory signal in the experimentally measured current. The 
time series were not long enough to capture many periods 
of the oscillation, if it was present. The narrow window 
of voltages over which the oscillatory regime would exist 
was not sampled sufficiently. New experiments, optimized 
for higher bandwidth and longer data acquisition of the 
current- voltage characteristics close to onset would be re- 
quired to detect the oscillatory regime, supposing that it 
occurs as predicted by the simulation. 

The laboratory measurements unambiguously show a 
sequence of large subcritical bifurcations as a function of 
1Z, for fixed IZe. As IZe is increased, these bifurcations 
evolve and proliferate in a regular way [14]. Based on 



measurements of the current, and some crude flow visu- 
alization, it was supposed that these were jump transi- 
tions between homogeneous azimuthal Fourier modes of 
the form m — > m ± 1. The size of the changes in the cur- 
rent were consistent with this supposition. Our numerical 
simulation of the same situation in which the flow pat- 
terns are revealed tells a different, more interesting story, 
however. 

As Fig. 4a and b show, the large subcritical transition 
shown in Fig. 2c takes the traveling, oscillatory vortex 
pattern to a strongly localized state consisting of just two 
isolated vortices. These vortices travel azimuthally in the 
direction of the shear, but are otherwise steady. Fig. 4b 
shows the perturbed electric potential associated with the 
localized vortices. The fluid flow within the vortices re- 
sembles that shown in Fig. 5d, which shows the perturbed 
stream function of a similar two vortex state under a larger 
shear. Fig. 4c shows the evolution of the mode amplitudes 
during the transition to the localized state, ra = 2 and 
several harmonics combine to make up the localized vor- 
tices, indicating the broadband nature of the flow, which 
clearly requires strong mode coupling to remain coherent. 

Upon reversing direction by decreasing 1Z in the range 
of 342 < 1Z < 388, the convection exhibits hysteresis as 
shown in Fig. 2c. Time series data of the kinetic energy for 
decreasing 1Z shows over-damped oscillations. While de- 
creasing 1Z from 500 down to 340, the originally dominant 
Fourier modes at ra = 2 and ra = 4 decay to essentially 
zero and m = 6 decays to a steady value. On decreas- 
ing 7£, the Nu — 1 curve shown in Fig. 2c reaches what 
is presumably a saddle- node bifurcation endpoint, where 
it rejoins the steady, non-oscillatory, m = 6 convection 
regime just above the onset at 7Z C . 

Fig. 5 shows a simulation which captures two successive 
subcritical bifurcations under a larger shear. In this case, 
the lower bifurcation results in a state consisting of a sin- 
gle traveling vortex. The subsequent transition leads to a 
two vortex state, as before. Several saddle nodes and new 
oscillatory localized states are seen. These successive tran- 
sitions in the simulation resemble the tree of bifurcations 
observed experimentally [14]. 

Thus, we find numerically that the mode structure 
changes near the strongly subcritical secondary bifurca- 
tions are much more complex than was apparent from the 
IV data alone. In particular, the modes may be oscilla- 
tory, and mode transitions can span more than ra — » rail. 
The traveling localized patterns above the bifurcation is 
not well characterized as having a single mode, and dur- 
ing the transition numerous coupled modes become active. 
The localized states also carry orders of magnitude more 
current than the traveling and oscillatory states below the 
bifurcation — so much so that the large subcritical jump 
was misidentified as the primary bifurcation in previous 
experimental studies [13,14]. 

Localized states in other systems [1] have sometimes 
been interpreted to be a consequence of heteroclinic cy- 
cles in the spatial coordinate [22] connecting a uniform 
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Fig. 4: (Color online) The spatial patterns of the electric po- 
tential perturbation, (i.e. the total potential minus the con- 
duction profile) at K = 388.0 and lie = 0.124, as in Fig. 2c. 
(a) Initially, there are six oscillatory traveling vortex pairs, (b) 
After 30r g , there are two non-oscillatory traveling vortices, (c) 
Time series of the Fourier amplitudes of the perturbed electric 
potential at midradius, showing the mode evolution during the 
transition. 

background state to a spatially modulated state and back 
again. Such states typically are stationary and occur in 
ID systems with a reflection symmetry. Here, such a sym- 
metry is inherently strongly broken by the handedness of 
the applied shear. In contrast, we find localized single vor- 
tices with only one sign of vorticity; they rotate internally 
in the same sense as the inner electrode, and travel az- 
imuthally in the same direction. The traveling aspect of 
the pattern may not be crucial, however, because it could 
be transformed away by passing to an appropriate rotating 
reference frame in which the localized state is stationary. 
The dynamics in the 2D annulus is in fact unaffected by 
overall rotations (see Ref. [12] for a detailed explanation 
of this rather counterintuitive fact). 

In the heteroclinic cycle model, the localized state can 
be thought of as two pinned fronts connecting states near 
the Maxwell point in a region of bistability between the 
two states. The Maxwell point is the parameter value at 
which the "energy", or Lyapounov functional, of the two 
states is equal. Near this point, the so-call snaking re- 
gion contains localized solutions with various numbers of 
internal spatial modulations. This picture has been ap- 
plied to the case of localized "convectons" in binary fluid 
convection for example [6,7], where the concentration and 
temperature fields are respectively stabilizing and desta- 
bilizing, somewhat analogous to shear and applied electric 
potential in our case. It is far from clear, however, whether 
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Fig. 5: (Color online) Numerical results for the dimensionless 
convective current, Nu — 1, at a = 0.56, V — 75.8 and He = 
0.231. In comparison with Fig. 2c for which He — 0.124, this 
larger shear produces more bifurcations over a similar range 
of the dimensionless control parameter. The spatial patterns 
of the perturbed stream function shown below are marked by 
letters in the upper plot, (a) Oscillatory convection ^XlZ — 600 
at t = 30 r q . (b) Undulating convection at a mode-changing 
bifurcation at 1Z — 650 at t — 30 r q . (c) Oscillatory and 
localized convection, consisting of a single vortex at 1Z — 800 
at t = 20 r q . (d) Non-oscillatory, traveling localized convection 
at n = 900 at t = 20 r q . 

this snaking picture can be relevant to the present case. 
We have not yet observed localized states consisting of 
multiple counter-rotating vortices, only isolated vortices 
of a single handedness surrounded by regions of very weak 
reverse vorticity. Also, the range of existence of the lo- 
calized states appears to span a range of parameters too 
wide to be confined to a narrow snaking region. They 
emerge not from a uniform state, but rather from oscilla- 
tory low-amplitude precursor patterns. Of course, neither 
the experiment, nor the time-stepping simulation, give us 
full access to all the solutions, stable and unstable, that 
may exist in the vicinity of the bifurcation. It is also 
worth noting that while the simulation reproduces the ex- 
perimentally measured currents quite well, the localized 
states have not been observed directly using the limited 
flow visualization available. 

For larger 1Ze and 7£, we observe the onset of chaotic 
flow, which followed the Ruelle-Takens-Newhouse scenario 
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Fig. 6: (Color online) Projections of the system attractor 
of Nu(t) — 1 using a time- delayed reconstruction method in 
sheared convection with a = 0.47, V — 16.3 and a large 
He = 0.8. (a) 2-frequency quasiperiodic flow (2-torus) at 
H — 639.3. (b) 3-frequency quasiperiodic flow (3-torus) at 
H = 673.4. These are spatially localized states which exhibit 
multifrequency oscillations in the time domain as the driving 
force is increased toward to a chaotic state. 

for the route to chaos [17]. Fig. 6 shows the time domain 
of Nu — 1 for a = 0.47, V = 16.3, and the large shear 
He = 0.8. With increasing 7£, we first find steady convec- 
tion, then a periodic state (a period-1 limit cycle), as in 
Fig. 3, followed by a two-frequency quasiperiodic flow (a 
2-torus), a three-frequency quasi-periodic flow (a 3-torus), 
and finally chaos. Fig. 6 illustrates the attractors of those 
states using a time-delay reconstruction [23] from time- 
series of Nu — 1. We first find a 1-frequency periodic mo- 
tion at H = 485.6 with a basic frequency of f\ ~ 0.60 r" 1 . 
Subsequently, the basic frequency of the periodic state 
changes slightly in the range 0.40 — 0.50 r" 1 as H in- 
creases. The 2-torus motion at H = 639.3 has fundamen- 
tal frequencies fi = 0.70 r" 1 and ji = 0.80 r" 1 . In space, 
the flow is localized to a single traveling vortex like that 
shown in Fig. 5c. The 3-torus state at H — 673.4 is also 
localized and has the basic frequencies / = 0.31, 0.69, and 
1.31, in units of r" 1 . Finally, at higher H = 1190, chaotic 
motion is observed, yet the traveling state remains local- 
ized in space. The complex time dependence comes from 
small blobs of charge which are chaotically emitted from 
the electrodes near the vortex. 

In conclusion, direct numerical simulation of sheared 
electroconvection provides an interesting new window on 
localized states, bifurcations and chaos in 2D fluid flows. It 
nicely complements previous experimental studies, which 
were mainly limited to current measurements, by allowing 
visualization of the basic fields. The simulations revealed 
an unexpected low amplitude state just above the onset of 
convection, and oscillatory states which have not yet been 
systematically investigated in experiment. A strongly sub- 
critical bifurcation makes transitions between these oscil- 
latory states to novel new states with localized traveling 
vortices. These carry much more current. Their localiza- 
tion is presumably due to some sort of mode interaction, 
but the exact nature of their stability is unknown. New, 
more sensitive experiments with better visualization will 
be necessary in order to study these states experimen- 
tally. Under stronger electrical forcing and higher rates of 
shear, we also observed a transition to localized chaotic 



convection via a Ruelle-Takens-Newhouse scenario. The 
relatively simple geometry of this system makes it an ideal 
place to explore ideas about localized, chaotic and turbu- 
lent states both theoretically and experimentally. 
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